Strong friction limit in quantum mechanics: the Quantum Smoluchowski Equation 
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For a quantum system coupled to a heat bath environment the strong friction limit is studied 
starting from the exact path integral formulation. Generalizing the classical Smoluchowski limit to 
low temperatures a time evolution equation for the position distribution is derived and the strong 
role of quantum fluctuations in this limit is revealed. 
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Quantum systems coupled to a heat bath environment 
can be found almost everywhere in physics and chem- 
istry H]. Transport processes in Josephson junctions [|| 
or electron transfer reactions in large molecules H are 
typical examples. What one aims to describe here is the 
effective dynamics of the relevant system degrees of free- 
dom, i.e. the reduced dynamics. While the corresponding 
classical theory is well established and based on Fokker- 
Planck equations, the formulation of dissipation in quan- 
tum mechanics is more complicated. In general a simple 
time evolution equation for the reduced density matrix 
does not exist [11 ; a formally exact expression for the re- 
duced dynamics in terms of path integrals is available, 
but the path integral expression in many cases cannot be 
evaluated even numerically. In the last years efforts have 
focused on the weak coupling regime where a description 
in terms of Master equations is possible [Q . The opposite 
limit of strong coupling has been left basically untouched 
1^. Here we study this limit. For the first time we show 
that for large friction the exact dynamics of a dissipative 
quantum system can be cast into a time evolution equa- 
tion for the position distribution, the so-called quantum 
Smoluchowski equation derived below in Eq. (13). 

In classical physics the Smoluchowski limit is well- 
known in all areas |g]. For large friction the Fokker- 
Planck equation in phase space reduces to a time evolu- 
tion equation in position space, the Smoluchowski equa- 
tion (SE). There, the basic condition is a time scale 
separation between relaxation of momentum and posi- 
tion which allows for an adiabatic elimination of the for- 
mer degree of freedom. Corrections to the Smoluchowski 
equation turn out to be algebraically small in the friction 
constant. Applications are countless; recent examples in- 
clude transport in systems with fluctuating barriers [Q 
or ratchets g as found in tunnel diodes Q and com- 
plex macromolecules [B|, and decay in periodically driven 
metastable potentials |0]. The quantum Smoluchowski 
equation should be important for the description of sim- 
ilar processes at low temperatures. 

Quantum dissipation-The standard way to de- 
scribe dissipation in quantum mechanics is based on 
system-freservoir-models with a total Hamiltonian H = 
Hs + Hr + Hj m. The reservoir Hr (heat bath) con- 



sists of a quasi-continuum of harmonic oscillators which 
are coupled bilinearly with the system Hs via the in- 
teraction Hj. Classically, this model leads back to a 
generalized Langevin equation for the reduced system. 
The quantum dynamics for the reduced density matrix 
follows from p{t) = TrR{&c£>{-im/K)W{fi) exp{iHt/h)} 
with the correlated initial state W^(0). In the position 
representation the path integral approach allows for an 
exact elimination of the bath degrees of freedom. Now, 
as in the classical Smoluchowski limit we focus on the 
position probability distribution P{q,t) = p{q,q,t) and 
obtain the exact result 0,|l^ 



P{<lf:t)^ / dqidq[Jt{qf,qi,q[)K{qi,q[). 



(1) 



Here, the propagating function Jt is a three-fold path 
integral (two in real, one in imaginary time) over the 
system degrees of freedom only. The real time paths q{s) 
and q'{s) run in time t from qi and q[ to the fixed end- 
point g/, while the imaginary time paths ^((t) connect qi 
with q[ in the interval Ti[3 (/3 = l/fceT'). The contribu- 
tion of each path is weighted by exp(iE[g, g,q']/?i) with 
an effective action E ||,|l^ not written down here explic- 
itly. The real time paths describe the dynamics of the 
system and the imaginary time paths specify the initial 
probability distribution 



P{q^,0) = [pfM.<l'^)K<l^.<l'^)] 



(2) 



where Pp{q,q') ~ (g|TrR exp(— /3-ff)|(7') is the reduced 
equilibrium density matrix while the preparation func- 
tion A describes the deviations from equilibrium. This 
way we avoid a factorizing initial state used in ordinary 
Feynman- Vernon theory J l|, which does not apply for 
strong friction. In Eq. (|l) the influence of the bath is 
completely determined by the spectral density /(w) of the 
bath oscillators. Effectively, the reduced system gains an 
additional interaction contribution (influence functional) 
which is non-local in time and strongly depends on fric- 
tion and temperature. Since the influence functional cou- 
ples real and imaginary time paths, an exact solution to 
Eq. (|^) is accessible only in certain cases. For weak damp- 
ing the Born-Markov approximation can be invoked lead- 
ing to the known Master equations 0]. Here, we explore 



to what extent the dynamics in (|l|) can be well approxi- 
mated by a time evolution equation when the friction is 
large. For a particle of mass M we assume Ohmic damp- 
ing /(w) = 2M7 (damping constant 7) supplemented by 
a high frequency cut-off Wc > 7 according to spectral 
densities in real systems and to avoid the well-known ul- 
traviolet divergencies of a pure Ohmic model. Further, as 
in the classical SE we take sufficiently smooth potentials 
for granted in this study. 

Harmonic oscillator- For the quantum Brownian mo- 
tion in the harmonic potential V{q) = \Muj'^q^ the prop- 
agating function is known exactly |lp^. To get insight 
in the Smoluchowski simplification we particularly study 
the relaxation of {q{t)) and {p{t)) in the overdamped 
limit 7/cjo ^ 1 and for times t 3> I/7. One finds 
{q{t)) « {q{t)),^ + 5q{t) and (p(t)) « {q{t)),,/^ + M5q{t)- 
The classical part shows the well-known sluggish decay 
{<l{t))c\ ex; exp(— Li;Qt/7) and to leading order {p{t)) « 
MSq{t). Quantum effects show up in Sq{t) which de- 
pends on the initial correlations between system and 
bath and for temperatures T > decays as exp{~i^t) 
where 1/ = 2tt /Tifi is the Matsubara frequency. Hence, 
on the time scale where the position relaxes, i.e. "i/oj^, 
the momentum can be considered as already equilibrated, 
i.e. {p{t)) « 0, only if a separation of time scales is 
guaranteed, i.e. ^/ujq ^ ?i/3, 1/7 [cf. Fig. 1]. This lat- 
ter condition defines the Smoluchowski range also for 
temperatures where quantum fluctuations are important. 
Now, neglecting terms of order exp(— 7^) and exp(— i^i) 
or smaller the time evolution for P{q,t) follows from 



dtPiq,t) = iu;l/j)d, [q+{q')d,] P{q,t) 

where the variance reads 

\ 2 1 



(3) 



(4) 



Inspection of the sum in Eq. (|4|) reveals that the Smolu- 
chowski range ^ /ujq ^ ?i/3, 1/7 comprises two particular 
subsets: one is the known classical range, 7 ^ i^, the 
other is a quantum mechanical region with 7 3> i^, un- 
studied so far (Fig. 1). While in both ranges to lead- 
ing order {q^) « {q^)c\, differences appear in the respec- 
tive corrections. Classically, these are of order fi (i/M 
and negligible against dynamical corrections of order 
l/M/37^. In contrast, in the quantum range to first order 



(?') - (<z')ci ~ {h/TTMj) log(?i/37/2^) = A 



(5) 



with A > l/M/37^. We conclude that for hj :^ ksT the 
dynamics of the quantum oscillator follows in leading or- 
der from the classical SE with quantum fluctuations as 
the dominating corrections. Since these are no longer al- 
gebraically small, they are of substantial relevance (see 
below). We note in passing that in the quantum range 



one flnds {p^} « (MH'y /n)[l + log(tJc/7)]- 
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FIG. 1 Smoluchowski range 7/1^0 ^ ^/^i I/7 (shaded). The 
classical range (7 <C i/) is simple shaded, the quantum range 
(7 2> 1^) double shaded. 

Quantum Smoluchowski egMation-Qualitatively, to 
leading order we expect the position diagonal element 
P{qf,t) to obey the classical SE even in anharmonic 
potentials and throughout the quantum Smoluchowski 
range (QSR) 



7/wo > np, 1/7 and ?i7 > fceT, 



(6) 



where ujq is interpreted as the ground state frequency. 
To quantify this statement we analyze the path integrals 
involved in Eq. (|^) in more detail. For that purpose, it 
is convenient to introduce difference and sum real-time 
paths X = q — q' and r = {q + q')/2, respectively. Then, it 
turns out that the action S][g, x, r] has friction dependent 
parts in imaginary and in real time |12| ] which suppress 
non-diagonal matrix elements. Accordingly, to leading 
order we assume \q{a) — q{0)\ and \x{s)\ to be sufficiently 
small and later verify this ansatz self-consistently. 

First, by expanding E[g, x,r] up to second order in 
the z-paths, we arrive at a solvable Gaussian x-path in- 
tegral where x and r paths are coupled bilinearly. We 
put x{s) — Xi + Sx{s) with fluctuations Sx{s) obeying 
6x{0) = Sx{t) = and to leading order obtain 



q^^^iisiq:=^,r]-T,lqfi,o])/h _ S{x,)e^'^^''^^^^'^'^''^ 



with the action 



S[r] = f ds[M-fr + V'{r)] 
Jo 



(7) 



(8) 



The path integral over the imaginary time path q with 
q{0) = q{Ti(i) — ri is easily evaluated. It yields a contri- 
bution proportional to exp[—/3I^(ri)]. Eventually, Eq. (|]) 
for P{qf,t) — p{qf,qf,t) reduces to 



P{qf,t)= Jdq,Giqf,t,q,)P{q,,0). 



(9) 



Here, the propagator stems from the remaining r-path 
integral and thus, is a sum over all paths r(s) with 
r(0) = qi,r{t) = qf where each contribution is weighted 
with exp{~S[r]/AM^kBT). However, this is exactly 



the path integral representation of the classical Smolu- 
chowski dynamics [|6[. Hence, to leading order the quan- 
tum SE looks like the classical one also for anharmonic 
potentials: dtP = {l/M^)dqLc\P with 



ici = V' + k^Td„ 



(10) 



where V {q) = dV{q)/dq. To justify our initial assump- 
tion we estimate that [q{a) — ^(O)]^ is at most of order A 
[see Eq. (P)] and x{s)^ at most of order H^P/Mjt. 

Of course, the most interesting issue is to specify 
contributions from quantum fluctuations to Eq. (Rw. 
However, before we do so let us first look at the sim- 
pler static problem and determine the equilibrium dis- 
tribution. In the relevant parameter range the path 
integral for Pii{q) = P/siq^q) is calculated by putting 
q{<j) = qma{<7) +y{a). The minimal action path qma with 
9ma(0) = qma^T-P) = 9 and its Euclidian action 5[g,na] are 
evaluated up terms linear in A. The result is S{q)/h = 
(3[V{q) — X(3V' {q)"^ /2]. The path integral over y is approx- 
imated by a Gaussian integral since (gma — y)^^O(A) 
and is solved by expanding y{a) = A X) yn.6^P(*^n^) 
with the boundary condition Ey„ = 0. As an important 
result wc obtain the unnormalized equilibrium position 
distribution in the QSR 

Pf3{q) = Pf{q) exp {a/3 [pV'{qf/2 - W"{q)/2\ } (11) 

with Pfjiq) = exp[— /3y(q)]. Next order corrections are 
much smaller, namely, of order {'h(3V" /M^Y . Amaz- 
ingly, for low temperatures P^ will be squeezed by fric- 
tion without limit. 

Let us now turn to the dynamical aspect of the QSR 
and evaluate quantum fluctuations to Eq. (IQ). For this 
purpose we write dtP = {l/M'-f)dqLq^P with Lqm = 
Lci + XSL where Ld is given in ( fiX\ j and 6L = SV' + SDdq 
with appropriate corrections to the potential and diffu- 
sion terms. To determine SV and SD requires to solve 
the full dynamical problem Eq. (mj perturbatively in A. 
Further, we need only to study the broad time range 
h(3, 1/7 <^ t < j/luq. Once we have SL in this time 
interval, we have it for all times as the time scale sep- 
aration renders the operator L^un time independent on 
the coarse grained time scale (^. The general strat- 
egy is as above: evaluate the minimal action paths to 
E[g, a:,r] and the corresponding fluctuations including 
corrections of order A. Again, a small x expansion of 
E allows us to neglect higher than bilinear couplings be- 
tween X and r-paths. Within the relevant time range 
it is also consistent to assume the difference \r{s) — r^ 
to be sufficiently small so that we may write y[r(s)] w 
Vin) + V'inMs) - r,] + V"{ri)/2[ris) - n]\ 

The solution to the classical a;-path is then simply 
Xc\{s) = Xieyi\i[—sV"{ri)/M"f]. Since the r-path is cou- 
pled to both the x and the q path, the expression for rd is 
rather lengthy and not given here. The same holds true 



for the imaginary time path, which can be put in the form 
qijj) = ri + Xqi [a] x] where qi contains also the coupling 
to the x-path. As for the equilibrium case, the fluctuation 
path integrals can be evaluated in a Gaussian approxi- 
mation: The contribution from the imaginary time path 
reads as in Eq. (|ll|), while the real time path integrals 
provide a factor of order 7. As a result we gain the prop- 
agating function in the QSR which determines P{q, t) via 
ordinary integrations over the initial coordinates Xi and 
Ti, see Eq. (|l|). Now, further progress can be made by 
introducing the scaled coordinate ki = Xi^. Then, as 
a function of ki the propagating function is basically a 
Gaussian with a width of order M^P/t, so that for suf- 
ficiently smooth preparations in x^-space, i.e. sufficiently 
localized preparations in momentum space, it is consis- 
tent to set A{ki/j, Ti) « A(0, r^). This allows us to carry 
out the ki integration. The time evolution of P{q, t) in 
the QSR and within the relevant time range is therefore 
of the form (H) with the propagator 



G{q,y,t) 



x/'ina{y,t) 



AQ~r{v,t)]''/4a{v,t) 



(12) 



Here, the width is ct = cro + Acti with ctq = {'t/Mj(3)[l — 
V"t/Mj] and cti = -[1-V"t/Mj], and the local motion 
reads r ^rQ + Ari with tq = y - [V't/M-/]{1 - V"t/2M-f) 
and ri = yf3V" <7i. 

The propagator (|l^) is now used to determine 6V' 
and SD in iqm. Starting from Eq. (|9|) one applies 
dt — {l/M^)dqLq„^ to G{q,y,t) and then expands y de- 
pendent terms in the integral up to second order around 
y — q. After tedious algebra we find 6V' — V'" /2 and 
SD = f3V". This leads us to the central result of this pa- 
per, namely, the dynamical equation of an overdamped 
quantum system, the so-called quantum Smoluchowski 
equation (QSE), dtP = {l/Mj)dgLq,aP with 



iqm = V' + W" 12 + dq fcBr[l + XfiV" 



(13) 



Two observations can be made here: First, quantum ef- 
fects give rise to an effective potential V^s = y + W" /2 
and an effective diffusion term Doff = k^T[l + \f3V"]. 
Second, seen as a continuity equation it turns out that the 
equilibrium distribution Eq. ( pT| ) indeed fulfills iqmP/3 = 
0. Corrections to the QSE are of order l3{V"hf /M^^'^ 
and can thus savely be disregarded. 

Now, from the QSE the quantum analog to the classi- 
cal Langevin equation is (in the Stratonovich sense [p3[) 



M-iq{t) + V'{q) + XV'" = v/T+^Sy" S,{t) (14) 

with Gaussian noise {£,{t)) = 0, {i{t)i{t')) = 
{2M^ / (3)S{t — t'). Hence, the quantum stochastic pro- 
cess in the QSR is equivalent to a classical process in an 
effective potential and with multiplicative noise. 

Applications-We illustrate the significance of quantum 
fluctuations in the QSR in two paradigmatic applications. 



Imagine a metastable barrier potential V{q) where the 
barrier height Vf, is the largest energy scale in the system, 
i.e. Vb ^ fceT'- Then, preparing the system initially in 
an equilibrium state restricted to the well region (around 
g = 0), for intermediate times (plateau range) the distri- 
bution becomes quasi-stationary P{q,t) -^ Pst{q)- This 
state corresponds to a constant flux across the barrier 
(around q = qb) Jst = (l/-^^7)-^qm^st determining the 
escape rate through T — Jst/N with the normalization 
N given by the well population. Now, adopting the clas- 
sical procedure we obtain from Eq. (p^) 



QSR 



\/V"{0)\V"{qb)\ ^_f3v, eA/3[y"(o)+iy"(g,)|] 



(15) 



Here, the second factor accounts for the quantum fluctu- 
ations and is not small at all (see Fig. 2). The point is 
that A dependent terms in the QSE enter exponentially 
and thus lead to a substantial rate increase compared to 
the classical rate. Already for moderate friction Tqsr 
agrees well with the exact rate (for Vb/fceT 3> 1) ||l4[] . 
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FIG. 2 Ratio of the exact rate with the classical (dashed) 
and the Smoluchowski (solid) rate vs. friction for ll/"(gi,)| = 
V"(0) = MujI 

In order to analyze V^ff and I?off in more detail we 
consider the dynamics of a Gaussian wave packet in a 
bistable potential V{q) = -Muj'^q'^ /2[1 - q'^/2ql]. The 
barrier top is located at g = and the wells at q — ±qa- 
Then, scaling times with j/loq and coordinates with 
a/S/Mwo we solve the time evolution of a wave packet 
initially localized near the barrier top numerically also 
for longer times. High precision results are presented 
in Fig. 3 for {q{t)) at two different temperatures. For 
moderate temperature the influence of V^ff prevails: the 
barrier height is diminished, i.e. Vcff^b = H ~ 3XMujq/2, 
thus increasing the population around g = so that 
{q{t)) < {q{t))c\. In contrast, at low temperatures diffu- 
sion effects dominate as the spreading of the wave packet 
near g = is retarded since the local effective tempera- 
ture is lowered Toff (0) < T, i.e. {q{t)) > (g(i))ci. 




Mean position vs. time of a Gaussian wave packet with 
{q{0)) ~ 0.5 moving in a bistable potential. Dashed (solid) 
lines are classical (quantum) results with A — 0.035. Lower 
(upper) curves are for uiohp = 1 (uiofiP = 6). 

To summarize we have derived the time evolution equa- 
tion for dissipative quantum systems in the strong fric- 
tion regime and have shown the substantial impact of 
quantum fluctuations. As the classical Smoluchowski 
limit is of importance for all overdamped systems at high 
temperatures, our results have potentially an equally 
wide range of applications for overdamped quantum sys- 
tems in low temperature physics and chemistry. 
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